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We investigate spin-one color superconductivity of a single quark flavor using the Ginzburg- 
Landau theory. First we examine the classic analysis of Bailin and Love and show that by restricting 
to the so-called inert states, it misses the true ground state in a part of the phase diagram. This 
suggests the use of the more general, noninert states in particular within three-flavor quark matter 
where the color neutrality constraint imposes stress on the spin-one pairing and may disfavor the 
symmetric color-spin-locked state. In the second part of the paper we show that, in analogy to some 
ferromagnetic materials, lack of space-inversion symmetry leads to a new term in the Ginzburg- 
Landau functional, which favors a spatially nonuniform long-range ordering with a spiral structure. 
In color superconductors, this new parity- violating term is a tiny effect of weak-interaction physics. 
The modified phase diagram is determined and the corresponding ground states for all the phases 
constructed. At the end, we estimate the coefficient of the new term in the free energy functional, 
and discuss its relevance for the phenomenology of dense quark matter. 

PACS numbers: 21.65.Qr, 74.20.De 

Keywords: Color superconductivity. Helical ordering, Ginzburg-Landau tlieory 



I. INTRODUCTION 

It has been known for a long time that certain mag- 
netic materials (e.g. MnSi and Fe^Coi-xSi) exhibit spa- 
tially nonuniform ordering with a long-wavelength he- 
lical structure. This phenomenon was first explained 
by Dzyaloshinsky [ij and Moriya fifty years ago. 
They pointed out that the lack of inversion symmetry 
of the underlying crystal lattice allows a new term in 
the Ginzburg-Landau free energy functional (hereafter 
referred to as the DM term), which makes the uniform 
ferromagnetic ground state unstable with respect to the 
formation of the helical spin density wave. A similar ef- 
fect also occurs in antiferromagnets Microscopically, 
the DM term results from the spin-orbit coupling. The 
theory of helical magnetism was developed to its present 
form two decades after the original discovery 0, Q , and 
is still a subject of intensive research in the condensed- 
matter-physics community @, 0, H] ■ 

In the present paper we show that this phenomenon 
has an analogy in a vastly different branch of physics, 
namely in strongly-interacting cold dense quark matter, 
which exhibits color superconductivity (see 3 for a re- 
cent review). The necessary ingredients for the helical 
ordering to take place are a vector (spin) order parameter 
and lack of space-inversion symmetry. We will therefore 
concentrate on spin-one color superconductivity. We will 
show that in this case, the parity-violating DM term is 
induced by weak interactions. 

Due to the exchange of gluons, the quark Cooper 
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pairs are formed predominantly in the spin-zero color- 
antitriplet channel. At very high baryon density where 
the quark masses can be neglected, the cold three-flavor 
quark matter is known to be in the so-called color-flavor- 
locked phase. However, at densities corresponding to the 
neutron star cores the large value of the strange quark 
mass as well as the charge neutrality constraint induce 
a mismatch of the Fermi levels of different quark flavors, 
and thus impose stress on the cross-flavor pairing. Other 
forms of pairing are then likely to occur. Depending on 
the size of the Fermi surface mismatch, quarks of two 
flavors and two colors may combine in the so-called 2SC 
phase. When the mismatch is too large, only quarks of 
the same flavor can pair and the spin-one pairing then 
remains the only possibility Q. 

Originally spin-one color superconductivity was sug- 
gested and studied as a mechanism for pairing of quarks 
of a single flavor or a single color left 

over from the 2SC pairing. The classiflcation and phys- 
ical properties of several spin-one color-superconducting 
phases were worked out in a series of papers by Schmitt et 
al. [l^, [H, [13, [H, EH ■ Possible impacts of spin-one color 
superconductivity on the phenomenolo gy of com pac t 
stars were studied in Refs. 0, [SI, [H [H, M [H, [11,113 ■ 
A different approach to spin-one color superconductiv- 
ity, based on the Schwinger-Dyson equations, was inves- 
tigated in [H]. 

The weak-coupling quantum chromo dynamics (QCD) 
calculations at asymptotically high density show that the 
ground state of the single-flavor quark matter is the color- 
spin-locked (CSL) phase. However, when spin-one color 
superconductivity is considered within three-flavor quark 
matter (e.g. as a complement to the primary 2SC pair- 
ing), the requirement of overall color neutrality may favor 
other patterns of spin-one pairing [2^ . We will therefore 
in the main body of this paper perform a phenomenologi- 
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cal analysis based on the Ginzburg-Landau (GL) theory, 
treating the coefficients in the free energy as unknown 
parameters. 

The plan of the paper is as follows. We start in Sec. 
im by investigating the possible forms that the order pa- 
rameter in a single-flavor spin-one color superconductor 
may take. Such a classification was already done by 
Schmitt [l9| for the sake of weak-coupling QCD calcu- 
lations. However, while in QCD the order parameter 
is provided by a nontrivial function of momentum near 
the Fermi surface — the gap function, in the GL descrip- 
tion we deal with a single (local) order parameter. This 
greatly simplifies the analysis, allowing us to carry out a 
complete classification of the possible symmetry-breaking 
patterns. 

In Sec. mil we review the GL theory for spin-one color 
superconductivity without the parity- violating DM term. 
This was already developed by Bailin and Love more than 
twenty years ago [13]. Nevertheless, they restricted their 
attention to the so-called inert states, proposing that one 
with the lowest free energy for the ground state. Here 
we derive a set of inequalities between the independent 
quartic terms in the GL functional which allow us to de- 
termine the unique absolute minimum of the free energy. 
We thus show that there is a sector in the phase diagram 
where the true ground state is actually noninert. 

In Sec. IIVI we finally introduce the parity-violating 
DM term into the GL free energy functional. Based on 
Ref. [§] we first work out the (slightly generalized) the- 
ory of the helical spin density wave in the ferromagnets. 
With all the necessary formalism ready, we then con- 
struct the corresponding helical states for the spin-one 
color-superconducting phases, and determine the modi- 
fied phase diagram. 

In Sec. |V] we demonstrate how the DM term arises 
from weak-interaction physics. We estimate the corre- 
sponding coefficient in the weak-coupling limit. In Sec. 
I VII we then discuss the possible relevance of this effect for 
the phenomenology of dense quark matter. In Sec. IVIII 
we summarize and conclude. 

It should be noted that while some calculations such as 
the derivation of the Ginzburg-Landau free energy in the 
high-density, weak-coupling approximation are standard 
and the details may thus certainly be omitted, the alge- 
braic analysis presented in Sections [TTl IHH and IIVI is par- 
ticular to spin-one color superconductivity. Even though 
all derivations are based on elementary linear algebra and 
require no lengthy computations, we choose to provide 
most details since they cannot be found elsewhere. For 
reader's convenience, we formulate some purely mathe- 
matical auxiliary material in the form of simple theorems 
and defer the full proofs to the appendices. 



II. CLASSIFICATION OF ORDER 
PARAMETERS 

In this section we will investigate the possible 
symmetry-breaking patterns in a single-flavor spin-one 
color superconductor. The order parameter transforms 
as an antitriplet under SU(3) color transformations, as a 
vector under spatial SO (3) rotations, and carries charge 
of the baryon number U(l) group. It can be represented 
by a complex 3x3 matrix, A^i, which transforms as (l9j 

A ^ UAR, (1) 

where U € SU(3) x U(l) = U(3)l and R G S0(3)r. The 
indices L and R denote the "left" and "right" symmetry 
groups, acting on the order parameter. This symmetry 
structure is similar to that of the superfluid Helium 3 
(30| . However, since the symmetry group of the spin-one 
color superconductors is larger than that of the superfiuid 
Helium [which has another SO (3) instead of the SU(3)], 
the classification will be somewhat simpler in the present 
case. 

In the following it will be helpful to consider also trans- 
formations from the "diagonal" subgroup, S0(3)v, 

A ^ R'^AR. (2) 

The classification of the possible inequivalent forms of 
the order parameter will be based on the following two 
claims which are proved in Appendix [XI 

Theorem 1 By a suitable symmetry transformation, the 
order parameter can always be brought in the form, 

(Ai ias -ia2\ 
-ia^ A2 iai (3) 
ia2 -iai A3 J 

(with real parameters Ai,ai), being a Hermitian, 
positive-semidefinite matrix. 

Theorem 2 Let the order parameter have the form ([3]) 
and U e U(3)l, R G S0(3)v. Then 

UR^AR = A (4) 

if and only if UA = A and R^ AR = A. 

Eq. ([3]) represents the simplest form to which the order 
parameter can in general be cast. As could have been ex- 
pected, it contains six independent parameters: A com- 
plex 3x3 matrix has altogether 18 real parameters, 12 
of which can be fixed by a transformation from the 12- 
parametric symmetry group, G = U(3)l x S0(3)r. We 
will classify all special forms of the order parameter which 
leave some continuous subgroup of G unbroken. The 
analysis is greatly simplified by Theorem [2] which en- 
sures that one can separately investigate invariance under 
left unitary transformations from U(3)l, and diagonal or- 
thogonal rotations from S0(3)v. [The transformed order 
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parameter ([T]) can always be written as on the left-hand 
side of Eq. ^ by the substitution U UR^ , which is 
just another matrix from U(3)l-] There is no nontrivial 
unbroken combination of transformations from the two 
groups. Therefore, we just need to classify the unbroken 
subgroups of U(3)l and S0(3)v- 

As follows from the polar decomposition, Theorem [3 
given in Appendix [XI U(3)l has a nontrivial unbroken 
subgroup if and only if the matrix A has zero modes. 
Specifically, the unbroken subgroup will be U(n)L where 
n is the number of zero modes of A. 

The possible unbroken subgroups of S0(3)v rnay be 
found by elementary geometry. The Hermitian order pa- 
rameter A is written as a sum of its real symmetric and 
imaginary antisymmetric parts, A = S + iA. These obvi- 
ously transform separately under the diagonal subgroup 
S0(3)v. Moreover, the antisymmetric part is parame- 
terized as Aij = eijkdk so that the three components 
transform as a vector, a. The (real) symmetric matrix 
S may be viewed as defining a quadratic surface with 
principal values Ai, A2, A3. Apparently, it possesses the 
same symmetry as A if and only if this quadratic sur- 
face is axially symmetric, with the axis given by a. This 
constrains S to be of the form 

Sij = aSij + PaiGj, 

or equivalently a linear combination of identity and the 
projector to the plane perpendicular to a, Vij = Sij — 
2 . Since S is actually diagonal, obviously at most one 
component of a may be nonzero (or (3 = 0) in order to 
preserve a continuous subgroup of S0(3)v. 

By combining the conditions for invariance under the 
left unitary and the diagonal orthogonal transformations, 
one arrives at the classification summarized in Fig. [1] 
The logic of the table is simple. First four rows display 
phases with nonzero a in order of increasing unbroken 
symmetry. The next two rows show the CSL and polar 
phases, whose order parameter is diagonal. Then come 
the phases Ni,N2 that completely break the S0(3)v. 
These cannot be easily cast in the form ([3]) and we thus 
display them in such a way as to manifest the number of 
zero modes. Finally, the axial and planar states, showed 
under the double Hne, are just special cases of the oblate 
and cylindrical ones, already included in the table. They 
are distinguished by an unbroken discrete symmetry, gen- 
erated by the permutation matrix 

/O 1 0\ 
P= 1 , 
\0 -1) 

that is, a rotation by tt about the line x = y, z = 0. We 
do not have the ambition to extend the above analysis to 
unbroken discrete symmetries. These two special cases 
are mentioned explicitly because they will later turn out 
to occupy a part of the phase diagram. 

Four of the indicated phases are "inert", i.e., their or- 
der parameter is fixed up to a symmetry transformation 



and an overall normalization. These are the A, CSL, pla- 
nar, and polar phases. Note that some of the states that 
are distinct according to Ref. [l^ , are classified as equiv- 
alent here. This is because (as already remarked above) 
we treat the matrix elements of A as pure numbers, not 
functions of momentum. 

A remark about the nomenclature is in order here. The 
names of the four inert phases are standard in literature 
on spin-one color superconductivity. The terms "oblate" 
[30I I and "e" [31] have been taken over from literature 
on superfluid Helium 3. The remaining four names are 
new. The "N" states are labeled by the degeneracy of 
the zero eigenvalue of A. In the "cylindrical" phase, the 
quadratic form S defines a cylinder with the vector a 
pointing along its axis. Finally, the "axial" state may be 
thought of as deformed CSL with just axial symmetry. 
As wc will see in Sec. IIV B 3[ this exactly happens to the 
CSL phase upon switching on the DM term [49l |. 



III. GINZBURG-LANDAU THEORY WITHOUT 
DM TERM 

After working out the classification of all inequivalent 
forms of the order parameter, we now investigate using 
the Ginzburg-Landau theory, which of the states is actu- 
ally energetically preferred. Our analysis will be similar 
to that of Baym and lida for spin-zero color supercon- 
ductivity f33|. So far, we do not include the DM term. 
In the following one should always keep in mind that the 
GL theory is strictly speaking only valid near the crit- 
ical temperature. For the sake of brevity, we often use 
the term "ground state" where we mean the state mini- 
mizing the free energy, i.e., thermodynamic equilibrium. 
It is only in Sec. IVII that we make some speculations 
concerning the physics far from the critical point, at low 
temperatures. 

It will sometimes be convenient to use a different no- 
tation for the order parameter, in particular to treat the 
matrix Aaj as a collection of three complex vectors, (jja, 
one for each anticolor a. The GL free energy density up 
to fourth order in A then reads 

+ b^l- $a+diA+d2B + d3C. (5) 

For the time being, the parameters ai, 02, b, di,d2, d^ are 
treated as free, constrained only by the requirement of 
boundedness of the free energy from below. There are 
three independent quartic G-invariant terms that we de- 
note as A, B, and C. Their explicit forms in both nota- 
tions for the order parameter are summarized in Fig. [21 
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FIG. 1: Classification of all order parameters that leave some continuous subgroup of G unbroken. The parameters Zi are 
complex, while all other parameters are real, in accordance with Eq. @. For convenience, we explicitly distinguish the SO(2) 
group of real rotations from the U(l) group of phase transformations, even though the two are actually isomorphic. The first 
eight rows represent all possible phases distinguished by the unbroken continuous symmetry. The last two rows — the axial and 
planar states — are merely special cases of the oblate and cylindrical ones which have an additional discrete symmetry. 
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FIG. 2: Independent quartic invariants and their expression 
using the matrix A as well as the complex vectors (f)a- 



A. Ground state 

Since the gradient part of the free energy is required to 
be bounded from below, the ground state is apparently 
a uniform field configuration that minimizes the static 
part of Eq. As usual, the "mass term" h changes 



sign at the critical temperature. In the following we will 
assume that it is negative (i.e., we are in the superfluid 
phase) so that the free energy has a nontrivial minimum. 
We will use the invariant A to measure the size of the 
condensate; actually, it is a squared norm of A in the 
sense that will be specified later. The values of the other 
invariants i3, C then measure the orientation of the order 
parameter in the color and spin space, and we will draw 
the phase diagram in the two-dimensional space of the 
parameters ^2,(^3. 

Trying to determine the ground state by a direct so- 
lution of the gap equation would be just hopeless. First, 
the gap equation is a coupled set of equations for the 
six independent parameters in A. Second, even if we 
somehow managed to solve it, we could at best show 
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that the solution is a local minimum of the free energy. 
Instead, we derive a set of inequalities between the in- 
variants A, B, C that allow us to uniquely determine the 
absolute minimum of the free energy. We again formu- 
late these inequalities as simple theorems whose proof is 
given in detail in Appendix [BI (It is understood that all 
the following claims about the order parameter hold up 
to a symmetry transformation.) 

Theorem 3 The invariants A, B satisfy the following in- 
equalities, 



1 



-A<B<A. 



(6) 



The first inequality is saturated (i.e., an equality holds) if 
and only if the order parameter is of the CSL type. The 
second inequality is saturated if and only if the matrix A 
has rank one. 

Theorem 4 The invariants B, C satisfy the following in- 
equalities. 



0<C<B. 



(7) 



The first inequality is saturated if and only if the order 
parameter is of the A type. The second inequality is sat- 
urated if and only if the order parameter is real. 

Theorem 5 The invariants A, B, C satisfy the following 
inequality. 



-A<B + C. 



(8) 



The inequality is saturated if and only if the order pa- 
rameter is of the oblate type with A2 — \J -f . 

Theorem 6 Let C < ^A. Then the invariants A,B,C 
satisfy the following inequality. 



Va<Vc + Vb~~c. 



(9) 



The inequality is saturated if and only if the order pa- 
rameter is of the e type. 

With this set of inequalities at hand, at is straight- 
forward to determine the phase diagram in the {d2,d^) 
plane. We do so by finding a lower bound on the free 
energy and showing that this bound is saturated by a 
particular type of the order parameter. In all cases, the 
free energy density can be, after fixing the orientation of 
the condensate, written as 



T = b\/A + dA, 

where d is an effective quartic coupling, specific for the 
given phase. The ground state condensate and free en- 
ergy are then given by 



2d 



1 ■•F Q 



Ad 



(10) 



> 



^2 + ^3 > 0, 1^2 > ^3. Using Theorems |4] and [5l we 

get 



d2B + dsC = ^{d2+d3){B+C) 

+ ^{d2-d3){B-~C) > l{d2+d3)A. 

To saturate the inequality, we should have simul- 
taneously B -\- C — ^A and B — C, i.e., B — ^A. 
By Theorem |3l this is only satisfied by the CSL 
state. The magnitude of the condensate and min- 
imum free energy density are given by ([10]) with 
dcsL = di + 

• (^2 + 0^3 < 0, ^3 < 0. In this case. Theorems [3] and 
H yield 

d2B + dsC > {d2 + dz)B > {d2 + d3)A. 

The order parameter that saturates this bound 
should be real and have rank one, which is precisely 
the polar phase, with condensate and free energy 
determined by dpoiar = di + d2 + d^. 

• (^2 < 0, da > 0. Now we estimate the free energy 
using Theorems [3] and [3] as 

d2B + dsC > d2A. 

The ground state is the A phase. The effective 
quartic coupling in this case reads dA = di -\- d2. 

• c?3 > (i2 > 0. This case is subtle; it is here that none 
of the inert phases provides the absolute minimum 
of the free energy. Let us assume that C < ^A. We 
then use TheoremlHland the Cauchy inequality (|B1|) 
with ui = \J d2{B -C),U2 = ^ (^2 + dz)C, v\ = 
1 / \fd2, V2 ^ 1/ \J d2 + ^3 , to obtain 



d2B + daC = (i2(S - C) + (^2 -I- d3)C = u\+ n\ > 

{uivi+u2V2)^ _ iVB~C + VCf ^ ^2(^2 + ^3) 

2d2 + d3 



> 



A. 



The inequality is saturated if and only if the order 
parameter is of the e type, and vectors u and v are 
collinear. This fixes the order parameter to be 



VA 



a ia 
-ia a 
Pj 



(11) 



2d2 



2d2 



The effective quartic coupling in this case is d^ = 
di + ^^4Pt^- One easily checks that the initial 

^ 2d2+d3 

assumption C < ^Ais fulfilled for considered values 

of ^2,^3- 



B. Phase diagram 

The calculation of the ground state for different rel- 
ative values of d2 , given in the previous subsection 
is straightforward, but may not be entirely transpar- 
ent. Therefore, we complement it here by an elegant 
and powerful geom etric picture, first developed by Kim 
and Frautschi [33l [33| to analyze complicated Higgs po- 
tentials in models of grand unification. It will not only 
confirm our previous conclusions about the ground state, 
but also illuminate the nature of the inert and noninert 
states and the phase transitions between them. 

The idea is as follows. The quartic part of the free en- 
ergy can be thought of as depending on the squared norm 
of the condensate. A, and two dimensionless quantities, 
A2 = B/A and A3 — C/A. These specify the orienta- 
tion of the condensate in the color and spin space. For a 
uniform field configuration, the free energy density thus 
becomes 

Tste.t^bVA+idi+d2X2+d3X3)A. (12) 

The inequalities derived above show that the quantities 
A2, A3 cannot acquire arbitrary values. Instead, their val- 
ues for all nonzero 3x3 matrices will span some domain 
in the (A2, A3) plane, which we will refer to as the target 
space, in order to distinguish it from the parameter space 
of ^2, c?3- The shape of the target space is a property of 
the algebra of 3 x 3 matrices (and the symmetry group 
G), and is independent of the couplings d2,d3. 

The absolute minimum of the free energy (|12p can now 
be found by a consecutive minimization with respect to 
the "angles" A2,A3, and then the "modulus" VA. But 
since the free energy ()12p is linear in A2,A3, the mini- 
mum will simply lie somewhere on the boundary of the 
target space. Which point of the boundary will realize 
the ground state, depends on the coefficients d2,dz- For 
fixed values of J-sta.t and ^/A, Eq. (fT2|) defines a straight 
line in the (A2, A3) plane. For a too small value of J^Jtat, 
this line will not intersect the target space, i.e., there is 
no state with the desired value of the free energy. As 
•?^3tat increases, the straight line will shift parallel until 
for some .Fmin, it will for the first time touch the target 
space. The point of touch will then define the ground 
state. 

In the case of the spin-one color superconductor, the 
explicit form of the target space is plotted in Fig. [31 We 
emphasize once again that it is the points at its boundary 
that will appear in the phase diagram as ground states 
for some particular combination of ^2,^3. The states 
that correspond to the corners of the target space there- 
fore play a distinguished role. These are the three inert 
states, A, CSL, and polar. Since they are inert, they are 
represented by a single point in the target space. A non- 
inert state of a particular type will occupy some nontriv- 
ial domain, depending on its number of free parameters. 

Let us now be more specific about the boundary of the 
target space. The second inequality in Eq. Q defines 
the edge connecting the "A" and "polar" corners, which 
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FIG. 3: Target space in the (A2,A3) plane (shaded). The 
various bounding curves are determined by Eqs. ([6l)-([9l). The 
curve connecting the points "A" , "X" , and "polar" is defined 
by VC -I- \/B — C = VA. For given d2, the structure of the 
ground state is determined by that point on the boundary of 
the target space which minimizes the expression d2X2 4- da A3, 
i.e., by the nearest point of the target space when looking at 
it (in the plane) in the direction of the vector (^2,(^3). 

thus involves all rank-one order parameters. The polar- 
CSL edge follows from the second inequality in Eq. ([7]) 
and consists of matrices A that can be made real by a 
symmetry transformation. The CSL-X edge comes from 
Eq. ([5]). It is occupied by matrices of the oblate type 
with A2 — -\/A| + (?. Finally, the curved segment X- 
A is a consequence of Theorem |6] and is realized by the 
matrices of the e type. This completes the picture of the 
target space and elucidates the significance of the various 
inequalities. 

The analysis carried out in Sec. IIII Al vields the phase 
diagram, displayed in Fig. [4] Note that all phase tran- 
sitions but the one between the A and e phases are first 
order. This is also easily seen using the Kim-Frautschi 
plot of the target space (see Fig. A straight segment 
of the boundary connecting two corners such as polar 
and CSL causes the ground state to change abruptly at 
an infinitesimal change of the slope of the straight line 
(fT2)) (i.e., the ratio of d2 and ^3). On the other hand, 
as we move through the e phase towards the transition 
line to the A phase, the order parameter runs along the 
curved segment X-A until it eventually continuously en- 
ters the A phase. Moreover, we can easily check explicitly 
using Eq. ifTTj) that as ^2 ^ 0-I-, the e state continuously 
goes to the A state. 

The Kim-Frautschi plot also tells us which types of or- 
der parameters may coexist right at the first-order phase 
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FIG. 4: Phase diagram of the spin-one color superconductor 
in the {d2, ds) plane. The solid and dashed thick lines denote 
first- and second-order phase transitions, respectively. The 
phase boundaries are defined by straight radial lines at angles 
j,^,TT,^ with respect to the d2 axis. The cross indicates 
the weak-coupling prediction, see Eq. (|?7)l . 



transition lines. If the segment of the boundary of the 
target space connecting the two competing phases were 
concave, we would find just these two states. How- 
ever, since all the border lines corresponding to first- 
order phase transitions are straight, a much wider class 
of states can actually coexist. Without going into de- 
tails we just note that special relations between c?2 and 
da which define the phase transition lines, may bring in 
additional degeneracy in terms of an enhanced symmetry 
of the free energy [35| . 

Note that the phase diagram in Fig. 2] does not depend 
on the coefficient di . The only way it affects the problem 
is indirectly, by the requirement of boundedness of the 
free energy from below. In other words, a particular value 
of di will determine a region in the {d2, da) plane which 
is physically allowed. 

Finally, a comparison with the phase diagram calcu- 
lated in Ref. [l^l [and plotted in the (^2, ^2 + d^) plane] 
shows that our results agree with the exception of the e 
region, which was missed in [l0| . The strategy used here, 
based on lower estimates of the free energy together with 
the conditions for their saturation, ensures that we have 
really found the absolute minimum of the free energy. 
Apart from the phase transition lines, it is unique up to 
a symmetry transformation. 



IV. GINZBURG-LANDAU THEORY WITH DM 
TERM 

We are now ready to analyze the GL theory for 
spin-one color superconductivity including the parity- 
violating DM term. However, since the construc- 
tion of the helix-ordered state for some of the color- 
superconducting phases is a bit involved, we prefer to 
illustrate the idea and develop the argument on a simple 
toy example. 

A. Toy model: Complex ferromagnet 

Following closely Ref. Q, we consider the GL free 
energy density functional of the form (soj 

T[M] = aiV.M^ • V,a1 + a2(V • M^)(V • M) 

+ bM'' ■ M + cNP ■ (V X M) -f d{M^ ■ Mf, (13) 

where the term proportional to c is the DM term. (Note 
that it is real up to a total derivative.) Without lack of 
generality we will assume that c > 0. 

First of all we would like to stress that since the DM 
term will make the ground state configuration nonuni- 
form, we are not a priori allowed to simply minimize the 
static part of the free energy to determine the magnitude 
of the order parameter M (which we shall in this section 
refer to as the magnetization) . Instead, we will rely on an 
estimate of the free energy, this time including the space 
dependence of the order parameter. Assuming the space 
has finite volume with periodic boundary conditions, 
we expand the magnetization in Fourier modes, 

a1(x) = ^ mke*'"''. 

k 

Using the integral form of the Cauchy inequality, we may 
estimate the quartic part of the free energy. 




The inequality is saturated if ■ M is uniform, i.e., 
the magnetization has the same magnitude in the whole 
space. 

Decomposing the Fourier mode mk into components 
parallel ("longitudinal") and perpendicular ("trans- 
verse") to the momentum k, TOk|| and TOk_L, and the 
transverse component further into its real and imaginary 
parts, mk_L = Uk + the DM contribution to the free 
energy becomes 

/ d^yi • (V X a1) = 2r2 ^ k • (uk X 4). 

k 
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Using elementary geometry, this is estimated as 

k • (uk X i7k) > -|k| |uk X Vk\ > 

> -|k| |uk| 141 > -|k|^!^l±M = _l|k| l^kj^. 

The chain of inequalities is saturated if the real and imag- 
inary parts of mk_L have the same size, are perpendicular 
to each other as well as to k, and together form a left- 
handed orthogonal system of vectors. (For c < 0, it would 
be right-handed.) Together with the gradient terms, the 
DM term can thus be minimized as follows, 

J d=^x[ai|V,Al|2 + a2|V • + cM'' ■ (V x M)] > 
>nY, [{ai + a2)|k|Vklll' + (ai|k|' - c|k|)|mk±|'] • 

k 

Stability of the system with respect to longitudinal fluc- 
tuations requires oi -I- a2 > 0. The longitudinal mode 
TOk|| therefore always increases the free energy. On the 
other hand, the DM term, being linear in momentum, can 
outweigh the gradient term and make nonuniform, trans- 
verse field configurations energetically favorable. The 
minimum free energy is achieved when only modes with 
|k| = |k|,nin = c/2ai are included. 

Putting all the pieces together, we obtain the minimum 
free energy density as 

^ J d^^J'[M]> (^b- M + dM^. (14) 

The state minimizing the free energy will be purely trans- 
verse so that the 02 term in Eq. (|13p actually does not 
play any role. To specify the form of the ground state 
more concretely, recall that apart from being composed 
solely of Fourier modes with |k| = |k|niin, it also ought 
to have a spatially uniform magnitude of magnetization. 
The most general state satisfying this condition has the 
form 

M ^ arhi.e'^ '' + /3m^e-'''"', (15) 

with fixed momentum k and real coefhcients a, (3. 

Two special cases deserve particular attention. First, 
if a or (3 is zero, the ground state is a single (complex) 
plane wave. Second, if we require the order parameter 
to be real (which is the case of the ferromagnet as well 
as several of the spin-one color-superconducting phases) , 
we arrive at 

M = rhi^e"^'^ + rnl.e~'-^'^ = 2[ukC0s(k-x)-UkSin(k-x)], 

i.e., a real standing wave. The magnetization evolves 
along a right-handed helix with the axis defined by the 
vector k, and the wavelength A = A-Kai/c. This concludes 
the argument and reveals the nature of the ground state 
induced by the DM term. 



Several remarks are in order here. First, for generic 
(nonzero) coefficients a, (3, the state (1151) breaks both ro- 
tational and translational invariance [51|. However, there 
is a combination of a rotation about k and a simultane- 
ous translation along k which remains unbroken. This 
leads to a peculiar, strongly anisotropic behavior of the 
Nambu-Goldstone mode of the broken symmetry 

Second, the wavelength of the helical state is propor- 
tional to 1 /c, hence the weaker is the DM term, the longer 
is the scale of the helical ordering. This will be particu- 
larly important in the later application to spin-one color 
superconductivity where the c-term comes from weak in- 
teractions and is thus expected to be tiny. 

Third, according to Eq. , the DM term effectively 
lowers the coefficient b. Therefore, it leads to a slight 
increase of the critical temperature, and at a fixed tem- 
perature to a slight increase of the magnitude of the mag- 
netization. However, its most notable consequence is the 
formation of the nonuniform field configuration. 

Fourth, as long as we restrict ourselves to terms up to 
fourth order in A in the free energy, the ground state 
can be rigorously proved to be composed of a single 
plane wave (and possibly the counterpropagating wave), 
although the gradient terms are minimized by Fourier 
modes with momenta lying on a sphere of radius |k|i„in- 

B. Spin-one color superconductor 

The generalization of the GL functional (fT3| to the case 
of the spin-one color superconductor is straightforward 
and unique, 

-f b$l ■ $a + C$1 ■ (V X $a) + diA + d2B + d^C . (16) 

As above, we Fourier-decompose the order parameter 
field, 

(?,(x) = ^^,ke'''-, (17) 

k 

and define A4 = X^k '/'Ik ' Vak- Using the same argu- 
ment as in the preceding subsection, the free energy is 
estimated as 

^ J d'j,J^[$a] >(^b- M + dM\ (18) 

where the effective quartic coupling d for the various 
spin-one color-superconducting phases was defined in 
Sec. nil Al The absolute minimum of the free energy 
is achieved for purely transverse configurations with uni- 
form summed magnitude (j>\-(f>a such that for all colors a, 
the real and imaginary parts of ipak have the same magni- 
tude and together with k form a left-handed orthogonal 
system of vectors, and only modes with |k| = |k|min = 
c/2ai are included. We shall now construct the helical 
states for the phases that appear in the phase diagram 
in Fig. m 
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1. Polar and A phases 

A single anticolor participates in the condensation. 
The ground state can therefore be constructed in com- 
plete analogy with Sec. IIV Al For the polar phase, the 
order parameter is real, and we thus find a real standing 
wave. Choosing the coordinate basis so that the helix 
points along the z axis, the order parameter takes the 
form 



*polar 



^ 0^ 


ycosfcz sinfcz Oy 



The condensate magnitude and free energy are given by 
a modification of Eq. (jlOp. 



2d '■ 



M 



(19) 



(20) 



mm 

where 

and as before, dpoiar = "i -i- U2 "3- 

In the case of the A phase, the order parameter can also 
be cast in such a form that only one anticolor condenses. 
(The matrix A is then of course no longer Hermitian.) 
The nonzero vector 0q is forced to be "maximally com- 
plex" in the sense that C = 0. This leads to the single 
plane wave solution. 



A 



mm ^ikz 




The condensate magnitude and free energy follow from 
Eq. (dH) with Ja = di + ^2- 



2. e phase 

In this phase, the order parameter can be cast in a 
form that just two anticolors condense, which amounts 
to adding another orthogonal real vector to the above 
A-phase order parameter. However, three real vectors 
cannot be simultaneously orthogonal to each other, and 
still orthogonal to the momentum k. This means that 
the lower bound (|18p cannot be reached, or, the gradient 
and static parts of the free energy cannot be separately 
minimized for an order parameter of the e type. 

To see what the ground state will look like, we resort to 
the Ginzburg-Landau equation following from Eq. ([Tl 



■ 2di(j)aU 



i>a - a2V(V • 0a) + 60 
l-$b)+'2d2M$a-cfl) 



f C(V X $a) 
2d3$l{$a-$b)^0. 



(21) 



We are not going to solve this equation directly. After 
all, even if we did, we could not prove anyway that our 
solution was the absolute minimum of the free energy. 
Nonetheless, Eq. (PT|) will provide us with the necessary 
insight to make a heuristic guess about the form of the 
ground state. 

The starting assumption is that the values of the order 
parameter (jja at different space points are connected by 
a symmetry transformation. This is reasonable for other- 
wise the static part of the free energy density would not 
be uniform; it is hard to imagine how such a configuration 
could be even a local minimum of the free energy. 

The DM term in the "equation of motion" ([2T|) forces 
the vector 0^ to rotate about the direction of k. Also, 
this term does not mix colors. Hence it is plausible that 
the transformation which connects the values of the order 
parameter at different points is a pure spatial rotation. 
Naturally, the axis of rotation is to be identified with 
the direction of the momentum k. We can then choose 
the basis in the color space in such a way that for one 
anticolor, the vector (j)a is perpendicular to k and rotates 
about it (transverse mode) , while the vector of the second 
anticolor is aligned with k and static (longitudinal mode). 
This is also in accord with our general discussion since a 
longitudinal mode with nonzero momentum would cost 
energy. 

With the above argument in mind, we write down the 
helical state of the e type with k along the z axis. 



^/2( 



V 




























(22) 



Here k,a,l3, and V Amin are treated as variational pa- 
rameters with the constraint 4a^ +0^ — 1 enforced by 
normalization. Minimizing the free energy within this 
class of variational states, we find 



ds) 



rf2(2 + + ^3(1 + ' 

q2 ^ ^2 - jdi 

^2(2 + 0+^3(1+0' 

and the momentum is, as above, given by fc = c/2ai. 
Note that for f = these expressions reduce to the pre- 
vious result, Eq. (fTTj) . The condensate magnitude and 
free energy read 



d-2+d-3 
2d2+d3 



^ mm 



2di 
62 1 



J d2+d3_ 

"2 2d2+d3 



2^ 



2d2+d3 



C2 di+d2+d3 
^ 2d2+d3 



di 



^ 2d2+d3 



While the denominators of the large fractions contain the 
expected effective quartic coupling, (4 = di + d2 24^+43 ' 
the numerators differ from Eq. . This is because the 
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separate minimization of the gradient and static parts 
of the free energy could not be achieved so that both 
the condensate and the condensation energy are actually 
smaller than Eq. would have predicted. 

Apart from the variational minimization of the free en- 
ergy with the ansatz ((22)) . we also checked explicitly that 
this state with parameters fixed by the above expressions 
indeed solves the GL equation (I^T]) . This provides a de- 
cent evidence that we have found the genuine ground 
state. 



3. CSL phase 

Here the situation is even more severe than for the e 
phase. In that case, the order parameter contains even 
after normalization a free parameter, say, the ratio a/ f3. 
The effect of the DM term in the free energy is then 
accounted for by a slight shift of this parameter. How- 
ever, the CSL order parameter is rigid, there are no free 
parameters beyond the overall norm to adjust. As a re- 
sult, the CSL state simply turns out to be incompatible 
with the DM term: It is no longer a minimum of the free 
energy. 

To see this, let us assume [as in the discussion below 
Eq. ([2T|l ] that the order parameter has everywhere the 
CSL form with a fixed magnitude. Consequently, the 
matrix A is unitary up to a real coordinate-independent 
factor. This dramatically simplifies the GL equation ([2T|) . 
which becomes linear and separated for the individual 
colors, 

-aiVVa - a2V(V • $a) X $a + 6off0a = 0, 



preferred. This argument leads us to the ansatz 



with &off = b + 2(icsLV.4. Being linear, this equation 
can be solved independently for the static and rotating 
modes. In particular the existence of a static mode re- 
quires bcff = 0. This fixes the condensate magnitude to 
its size without the DM term. The rotating mode has to 
fulfill the condition fli I kp =c|k|. Therefore, even though 
there is a solution with helical structure, its energy is the 
same as that of the uniform CSL condensate. 

In order to resolve this problem, it is important to real- 
ize that the anticipated existence of a helical structure in 
the ground state implies breaking of the rotational sym- 
metry at least to the group of rotations about the axis 
of the helix, k. It would be naive to expect the isotropic 
CSL state in such a situation. 

We can take the reduced axial symmetry into account 
and at the same time relax the rigidity of the order pa- 
rameter by considering the more general axial state (see 
Fig. [T]). This is natural: There is no reason why the static 
part of the order parameter, aligned with k, should have 
the same length as the transverse part, which is perpen- 
dicular to k and rotates about it. In fact, thanks to the 
DM term, we should expect the transverse part to be 



'axial 



min 



a cos kz a sin kz 
—as'mkz a cos kz 
(3i 



(23) 



with the normalization constraint 2a^ + 0^ = 1. The 
variational minimization of the free energy now results 
in 



1 



a 



f rfl+rf2+rf3 

3 + 2^ 



1-2^ 



3 + 2^ 



The condensate magnitude and free energy are given by 
h 3 + 2^ 



T ■ — 

mm ~ 



2 3di -f d2 + ' 
4 3di +d2+d3 



These results show that a is always larger than f3 so that 
the rotating part of the order parameter is indeed favored 
over the static one, as predicted. In comparison to the 
CSL state this physically means that one can gain energy 
by making the condensates of the three anticolors slightly 
imbalanced in favor of the two anticolors that form the 
helical structure. This has the amusing consequence that 
the ground state is no longer color-neutral. The color- 
density imbalance is a sheer weak-interaction effect. 



4- Planar phase 

The magnitude f3 of the static part of the axial order 
parameter decreases with the sum d2 + ds until at d2 + 
d^ — 2^di it goes to zero. The axial order parameter 
reduces to the planar one. Since in the phase diagram 
without the DM term. Fig. |4l the boundary between the 
CSL and polar phases occurs at ^2 + ^3 = 0, we may 
expect the planar phase to interpose itself between the 
two at nonzero ^. 

The planar order parameter is explicitly expressed as 



^planar 



-All 



cosfcz sinfcz 
— sin kz cos kz 




(24) 



and its size and free energy are given by Eq. (jl9p with 



-^planar 



2 



C. Phase diagram 



The phase diagram, modified by the parity-violating 
DM term, is displayed in Fig. [5l assuming that di > 0. 
[For di < 0, a large part of the (^2,^3) plane would be 
excluded by the requirement of the boundedness of the 
free energy so that the resulting picture would not be very 
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polar 




FIG. 5: Phase diagram of the spin-one color superconductor 
including the DM term. The solid and dashed thick lines 
denote first- and second-order phase transitions, respectively. 
The positions of the solid lines are the same as in Fig. |3]and 
independent of the DM term, whereas the offset of the dashed 
lines is proportional to the parameter ^ defined in Eq. (PU)) . 



interesting.] The e-A and axial-planar phase transitions 
are of the second order. They are both characterized by 
vanishing of the respective /3 parameter. The rest of the 
phase transitions are first order. 

For the axial and e states, the minimum free energy 
suggested by Eq. ([18]) cannot be achieved. As a result, 
these two phases are pushed away from the phase dia- 
gram by the DM term. The first-order transition lines 
were determined by a comparison of the free energies of 
the individual phases. It was thus checked that the free 
energy is continuous all over the phase diagram, which is 
a necessary condition for the conclusion that we have not 
missed any possible intermediate phase and determined 
the phase diagram correctly. 

Finally, note that while the complete GL free energy 
(HH) is SU(3) X S0(3) X U(l)-invariant, the DM term 
has in fact an enhanced SO (6) x SO (3) symmetry un- 
der which the pairs of states A-planar and e-axial are 
degenerate. This has in particular the consequence that 
the e-axial transition line in the phase diagram (the CSL 
state being a special case of the axial one) does not shift 
when the DM term is switched on. 



V. MICROSCOPIC DERIVATION OF DM 
TERM 

In this section, we provide a microscopic derivation of 
the Ginzburg-Landau free energy. Again, a thorough 
analysis of this problem was performed by Bailin and 
Love long time ago. Nevertheless, since the calculation 



at the level of generality taken up in Ref. [1(| is rather 
involved, we make a number of simplifying assumptions 
that allow us to arrive at a simple formula in an efficient 
manner. 



A. GL functional from NJL model 

First of all, we fix the Dirac structure of the spin-one 
gap matrix to be simply 7^ . Such an order parameter has 
positive parity and in the ultrarelativistic limit describes 
purely transverse pairing of fermions, which was shown in 
Ref. II9I to be energetically preferred to the longitudinal 
pairing '52|. 

Second, we for simplicity disregard the gluonic fluc- 
tuations which would otherwise make the superconduct- 
ing phase transition first order [3^, and derive the GL 
free energy in the framework of the Nambu-Jona-Lasinio 
(NJL) model [13] ■ This is decently justified by the fact 
that the coefficients of the GL free energy are to a large 
extent universal, the only dynamics-dependent quantity 
being the critical temperature jl3]. The fact that the 
NJL model does not capture correctly the effects of the 
soft chromomagnetic gluons and thus predicts wrong 
asymptotic behavior of the gap and critical temperature, 
should therefore not matter as long as the critical tem- 
perature is appropriately adjusted. 

Third, we neglect, as usual, the antiparticle degrees 
of freedom. As mentioned above, in the ultrarelativistic 
limit this automatically projects out the transverse part 
of the order parameter. Here we take the advantage of the 
fact that the Dirac structure 7^ can be achieved within a 
NJL- type model with a contact, momentum-independent 
interaction. Finally, we neglect the admixture of states 
with higher angular momentum [38,, ^] . 

Given the above assumptions, the mean-field approxi- 
mation to the NJL model amounts to a theory of nonin- 
teracting quasiquarks in the background of the (possibly 
slowly- varying) order parameter (jjai, with the imaginary- 
time propagator in the Nambu space, ^ = {tp,ip'^)'^ , 
given by 



_ / (iujn - ^k)7oA,| 



(i^n + Ck)7oA^ 



(25) 



where we used the usual notation for ek = 'v/l? + ttt-^ 
and the energy with respect to the Fermi level, = Ck ^ 



^. Also, A^ 



1± ^70(7 -k + m) 



are the standard 



positive/negative energy state projectors, and 



(26) 



where (Af )bc = -iiabc 



The GL free energy density at weak coupling and in 
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the ultrarelativistic limit then becomes 



2407r4^r2 
Stt' 



7t\<Pa\' 



7C(3)^Pf 
^ 4807r4T2 



(3y^ + 36-2C), (27) 



where Tc is the critical temperature for the spin-one 
pairing. We denoted t = ^ — 1 and for the sake of 
an easy check with literature, kept separate symbols for 
the Fermi momentum pp and the chemical potential fi, 
even though they actually coincide in the ultrarelativis- 
tic limit. Derivation of the GL functional ([27|) is routine. 
Yet we are not aware of any other independent calcula- 
tion for spin-one color superconductivity using the NJL 
model, and therefore provide some details in Appendix 

o 

A comparison with Eq. ([5]) now yields the phenomeno- 
logical coefficients 



oi = — 2a2 = 



1207r4^r2 ^ 



5 = 



37r2 



di= d2 = --da = 



7C(3)/ipF 
1607r4T2 ■ 



The ratio of d2 and ds is indicated in Fig. [J] by a cross. 
In agreement with the full QCD calculation we conclude 
that in the absence of the DM term, the weak-coupling 
limit favors the CSL pattern. 

In the presence of the DM term, the form of the ground 
state depends on the parameter 



, 907r6T2 
'7C(3)p|r 



As remarked at the end of Sec. IIV Al the DM term 
slightly increases the critical temperature. Eq. (fT9)) sug- 
gests that the corrected critical temperature is given by 
C = -l [53, i.e.. 



tr 



, 907rfar2 



Just below the (new) critical temperature, ^ is large and 
negative. With decreasing temperature, it goes through a 
singularity to large positive values and starts decreasing. 
(Of course, this "singularity" is completely artificial and 
comes just from the definition of ^.) In the tiny window 
below the critical temperature and above the tempera- 
ture at which ^ = 1/6, the ground state has therefore 
the planar structure, Eq. ([24]) . Only for ^ < 1/6 the 
system enters the distorted CSL phase — the axial phase, 
the ground state being as in Eq. (l23l) with 



1 + 4^ 
3 + 2^ 



1-6^ 
3 + 2f 



With further decreasing temperature, ^ drops to zero and 
the ground state relaxes to the CSL state. Needless to say 




FIG. 6: Pairing field self-energy with Z-boson insertion. 



that the temperature range in which this evolution occurs 
is extremely narrow, but it is nevertheless interesting to 
observe that the cooling of the system across the critical 
temperature actually consists of a fast sequence of two 
phase transitions, both of second order. 



B. DM term 

The last, and very important, missing ingredient in 
our analysis is the actual value of the DM coefficient c. 
We will argue here that the parity-violating DM term is 
naturally induced by the underlying weak interactions. 
Since the DM term is bilinear in the order parameter, we 
will seek weak corrections to the collective pairing mode 
propagator. 

Before we begin the calculation, we make a remark 
about the discrete symmetries of the DM term. As al- 
ready stressed, it breaks the parity. At the same time, 
it is invariant under charge conjugation (up to a total 
derivative). On the other hand, it is well known that 
weak interactions in a single-fermion-family world violate 
parity, but preserve the combined CP transformation. In 
the vacuum, the DM term would therefore be prohibited, 
at least as a consequence of weak interactions. However, 
charge conjugation is broken explicitly by the presence 
of the dense medium, and the DM term therefore arises 
from the interplay of weak interactions and many-body 
physics. 

Taking into account Gaussian fluctuations above the 
mean-field Cooper pair condensate, the pairing field 
propagator is given by a geometric series of the fermion 
bubble diagrams [i^]. The most straightforward weak 
correction is then the .Z-boson exchange, depicted in Fig. 
ini {W^ bosons cannot be exchanged since we consider 
just single-flavor quark matter.) Unfortunately, this does 
not work for the following reason. As mentioned above, 
the dominant spin-one pairing pattern in the ultrarela- 
tivistic limit is transverse in the relative momentum of 
the pair. Roughly speaking, it corresponds to pairing 
with total momentum and orbital momentum zero and 
total spin one, i.e., opposite helicity (and chirality). The 
quark interaction vertex with the Z boson can be param- 
eterized as 7'^(u -|- W75), where 



-— -^TT— sin^ 6'w - cos^ Ow ] , 
2sm20w \3 



2 sin 20W 
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q + k 



bj 



-q+p 

FIG. 7: Effective coupiing of tlie pairing mode to tlie Z boson. 
Tfie labels of the fermionic lines denote momentum flowing 
in the direction of the arrows. 




FIG. 8: Pairing field self-energy with Z-boson insertion. The 
hatched circles denote the effective Z-diquark vertex defined 
by Fig. m 



for u-type quarks, and 



2sin26'w I 3 

e 



sin^ 9w — cos^ 9w I , 



2 sin 29- 



w 

for d-type quarks. Here 6'w is the Weinberg angle and e is 
the electric charge unit. The coupling of the Z boson to 
two quarks of opposite chirality in the loop in Fig. [5]then 
produces the factor {u + v){u — v) — v?' — . However, 
to achieve parity breaking, one needs an amplitude odd 
in V. 

We therefore have to look for corrections of higher or- 
ders. Since each propagator of a weak intermediate bo- 
son is suppressed by the huge weak scale of the order of 
100 GeV, it is most likely that any diagram with the ex- 
change of more than one heavy vector boson will be much 
smaller than a graph with just one Z boson, but includ- 
ing fluctuations of the pairing field. Having in mind that 
the underlying interaction of the Z boson with the weak 
neutral quark current induces an effective coupling of Z 
to the pairing field 0a (see Fig. [7]), we thus anticipate 
that the leading contribution to the DM term will be 
given by the graph in Fig. [H The fact that the diagram 
involves propagation of an intermediate spin-one collec- 
tive mode unfortunately means that we will be able to 
provide only a rough, order-of-magnitude estimate of the 
DM coefficient c. 

We first focus on the effective Z-diquark vertex, 
r^^jj (fc,p), with the vector and axial- vector parts de- 
fined by 



r 



;_,^.(fc,p)=ziV;%^.(fc,p) + z;A^,^.(/c,p). 



For the sake of an analytical estimate of the DM term, 
we approximate these vertex functions by their Taylor 
expansion at zero external momentum, working again in 
the high-density approximation near the Fermi surface 
[41.] . For the vector part of the vertex we thus get 



2/i 



where Ap is an ultraviolet cutoff on the quark momen- 
tum; we will choose it at the order of the Fermi mo- 
mentum. On the contrary, the axial part of the vertex 
is finite; it yields the expected parity-violating structure 
and its Taylor series starts at the first order. Given the 
g^n structure of the Z-boson propagator, we just need to 
evaluate the axial vertex for /i = 0, 



7<(3)/iPF 
487r4T2 



The DM term is produced by the diagram in Fig. 
[8] where one of the effective vertices is vector-like and 
the other one axial-like. Neglecting the momentum- 
dependent part of the Z-boson propagator gives 



1 



Qck.di {k)g^,^ A^"^^^^^ {k - p, p)V^i,i,j {-k+p, k) 



+Valckik'P,PW,i,,,{-k+p,k) 



where we used the shorthand notation for a sum-integral. 



and Gai,bj denotes the collective mode propagator. In 
the normal phase, it is diagonal in the internal space, 
Qai,bj(k) = SabSijG{k). A Comparison with the DM term 
in Eq. lfT6|) then results in an expression for the coefficient 



c, 



7MwC(3)/i^PF 

367r6T2M| 



log- 



2AFeT-i 



irTr 



Gik). 



The last sum-integral is apparently quadratically di- 
vergent. Part of the divergence comes from the fact that 
we approximated the effective vertices in Fig. [8] with 
their low-momentum limits. We therefore make a rough 
estimate by replacing the Matsubara sums with a fre- 
quency integral and putting an ultraviolet cutoff, Ab , on 
the frequency-momentum integration. Taking into ac- 
count the appropriate volume measure and the fact that 
the coefficient of the leading, 0{p^), in the inverse prop- 
agator of the collective mode is typically of order 10~^ in 



the ultrarelativistic limit 421, we find 



7wwC(3)/x^Pf 



p« lO-^Al log 



2AFeT-i 



(28) 



The cutoff Ab should be well above the characteristic 
scale of the pair fluctuations in order not to suppress any 
physical contribution to the integral; for the moment, we 
treat is as a free parameter. 
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VI. PHENOMENOLOGICAL IMPLICATIONS 

With all coefficients of the GL functional we can now 
readily determine the parameters of the ground state. In 
order to be able to evaluate the condensation energy, we 
naively extrapolate the GL theory to zero temperature. 
For a concrete calculation we consider the CSL structure 
with d — di + ^2±d^ ^ in the ultrarelativistic limit where 
Pf = /i; the orders of magnitude will nevertheless be the 
same for all phases. 

First of all, the wavenumber of the helix in the ground 
state is given by 



«DM 



2ai 



To appreciate the robustness of the helical ordering, we 
next compare the condensation energies of the uniform 
condensation and of the helical ordering itself. The first 
reads 



1, 



Ad 



7C(3) 



where &o = h/t, while the latter is 



DM 



'4ai 9^6 \^M| 



[For simplicity, we use the lower bound from Eq. (jlSp 
which is not really saturated for the CSL state. How- 
ever, this only changes the result by a factor of order 
one.] Apparently, the energy gain from the formation of 
the helical structure can be many orders of magnitude 
smaller than the condensation energy of the supercon- 
ducting state. It is therefore natural to ask whether the 
long-range helical ordering will not be destroyed by ther- 
mal fluctuations of the order parameter. 

At low temperatures these will be dominated by the 
Nambu-Goldstone boson(s) of the spontaneously broken 
symmetry. The energy density deposited in the thermal 
fluctuations will therefore be that of the phonon gas with 
the phase velocity Wph = 1/ \/3 (in the ultrarelativistic 
BCS hmit), i.e., Stt^T'^ /90v^^^. The hehcal ordering will 
be destroyed when this becomes comparable to the con- 
densation energy, J-um- This leads to the characteristic 
temperature 



TT 



/j_ / uvp \ 

2 [AqJ 



1/2 



[We dropped the unimportant numerical prefactor 
(50/9\/3)i/4 which is close to one. 

For temperatures in the range Tdm ^ T < Tc, the 
system will look just like the uniform spin-one color su- 
perconductor. Only at temperatures lower than Tqm the 
helical ordering will appear. This change of behavior will, 
of course, not be associated with any phase transition. 



In order to assess the importance of the DM effect, we 
now make a specific numerical order-of-magnitude esti- 
mate. To that end, we need to know that the electroweak 
couplings are such that uv ^ 10~^, the Z-boson mass is 
of order Mz ~ 100 GcV. We also choose the typical 
value of the chemical potential to be /i = 400 MeV and 
set the fermionic cutoff equally, Ap = fi. There is some 
controversy in literature regarding the size of the critical 
temperature, or gap, in spin-one color superconductors. 
Let us be rather optimistic and assume that ~ 100 keV 
[Hi. Using Eq. ^ we thus get 



M 



MeV 



IQ- 



VMeV 



DM 



IQ-^Ab. 



The choice of the bosonic cutoff Ab is obviously crucial. 
Our rough guess is slightly above the critical temperature 
itself, Ab ^ 1 MeV. For this value the temperature scale 
at which the helical ordering takes place becomes 0.1 keV, 
and the helix wavelength is truly macroscopic — about a 
millimeter. The temperature scale suggests that the phe- 
nomenon could take place in extremely cold quark matter 
only (six orders of magnitude below the Fermi tempera- 
ture) such as in very old neutron stars. 

In any case, however, the helical ordering in the ground 
state results in peculiar low-energy properties of the sys- 
tem. As already remarked above, this nonuniform or- 
dering breaks both translational and rotational invari- 
ance, leaving unbroken just their special combination. 
The system will therefore exhibit anisotropic behavior. 
In helical ferromagnets, the Nambu-Goldstone boson as- 
sociated with the broken symmetry acts as a magnon 
when it propagates transversely to the helix axis, and as 
a phonon when it propagates along Q . The case of spin- 
one color superconductors will be similar, in particular 
we expect that the Nambu-Goldstone mode will have dif- 
ferent phase velocities in the transverse and longitudinal 
directions. This will in turn affect the thermodynamic 
properties of the system at very low temperatures. 



VII. SUMMARY AND DISCUSSION 

In this paper we have revised in detail the structure 
of the ground state of spin-one color superconductors, 
composed of a single quark flavor. In Sec. |TT] we pro- 
vided a complete classification of possible ground states 
distinguished by the unbroken continuous symmetry, as- 
suming a local constant order parameter. In Sec. IIIII 
we worked out the Ginzburg-Landau description of the 
spin-one color superconductor near the critical tempera- 
ture. We wrote down the most general GL free energy up 
to fourth power in the order parameter and determined 
its unique global minimum. We thus revealed that the 
noninert s state is favored in part of the parameter space. 

In Sections IIVI and |V] we then argued that the ground 
state of a spin-one color superconductor will actually be 
nonuniform as a result of a tiny parity-violating effect 
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which is due to the electroweak physics. Within the 
GL analysis, this effect is easily taken into account by 
adding a new term to the free energy. In analogy with 
some ferromagnetic materials, the ground state is found 
to exhibit helical ordering, typically with a wavelength 
much longer than the characteristic scale of the underly- 
ing many-fermion system. This leads in particular to an 
anisotropic behavior of the system at low temperatures 
where the helical ordering sets. 

In the course of our analysis, we made a number of 
simplifying assumptions that we wish to discuss to some 
extent now. First, we assumed that the local structure 
of the Cooper pair and hence also the classification of 
Sec. HI] does not change when the DM term is switched 
on. This is plausible when the wavelength of the helix is 
long enough, in particular much longer than the Cooper 
pair size, so that the spatial modulation of the condensate 
can be treated as a perturbation. In the realization of the 
phenomenon considered in this paper, this is guaranteed 
by the huge ratio of the electroweak vector boson mass 
to the pairing scale. 

Second, while we considered just a positive-parity spin- 
one condensate, electroweak interactions also induce a 
small admixture of a negative-parity condensate which 
makes the GL analysis more involved [l^. This conden- 
sate does not interfere with the DM mechanism to the 
first order in the electroweak effects, and we therefore 
neglect it here. Moreover, while this condensate would 
be extremely difficult to detect, the helical ordering of 
the ground state provides a macroscopic realization of 
parity violation. 

Third, there are other sources of parity violation than 
direct weak-interaction effects. For instance, condensa- 
tion of pseudoscalar mesons due to finite chemical poten- 
tial also breaks parity [45|. However, this is most likely 
to happen only in the color-fiavor-locked phase where all 
quarks are paired, and hence it is not relevant for spin- 
one color superconductivity. 

Fourth, the fact that the condensate is spatially 
nonuniform means that the quarks pair with a small, 
yet nonzero total momentum. In this sense, the heli- 
cal state is analogous to the Larkin-Ovchinnikov-Fulde- 
Ferrell state considered in crystalline color superconduc- 
tors [3, [1^. Different is, of course, the mechanism of 
the inhomogeneous pairing: In crystalline color super- 
conductors, it is driven by a mismatch of the Fermi sur- 
faces of the quarks to be paired, whereas in the present 
case the Fermi momenta are exactly equal. The nonzero 
momentum of the pair is induced by the DM term. The 
important difference between the various phases consid- 
ered here is that unlike the real axial, planar, and polar 
states, the single-plane- wave states A and e carry nonzero 
(charge as well as color) current. While in crystalline 
superconductors this is balanced by a backflow of the 
unpaired fermions so that the total net current in the 
ground state is zero the same issue in the A and e 
helical phases yet remains to be clarified. 

We observe that the axial and planar states which in 



a sense interpolate between the CSL and polar phases, 
show up as physical ground states in the phase diagram. 
This suggests that they should be seriously taken into 
consideration even in the absence of the DM term. In 
particular, in Ref. [29| the effect of color neutrality on 
spin-one color superconductivity as a complement to the 
primary 2SC pairing was studied. It was shown that the 
polar phase may be energetically preferred. Our results 
open the possibility of an even more favorable state of the 
axial type, which can compensate the color imbalance of 
the 2SC pairing and yet gain energy from condensation 
of all three anticolors. 

Finally, we emphasize that the helical ordering by the 
DM interaction is a fairly general phenomenon; the only 
two prerequisites are a vector order parameter and bro- 
ken parity. There are several issues that deserve further 
investigation. Besides the stress on the spin-one color 
superconductivity which stems from the pairing of the 
other quark fiavors, we would like to point out in par- 
ticular the peculiar properties of collective excitations. 
They will be anisotropic even in the CSL-like axial state 
and may thus in principle provide a clear manifestation 
of the nonuniform nature of the ground state. We are 
going to study these issues in our future work. 
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APPENDIX A: ORDER PARAMETER 

In this appendix we prove Theorems[l]and[2l The basic 
ingredient will be the polar decomposition, well known 
from linear algebra. 

Theorem 7 (Polar decomposition) Let M he a 

square complex matrix. There is a unitary matrix U and 
a positive-semidefinite Hermitian matrix H such that 

M = UH. (Al) 

The matrix H is unique, while the matrix U is unique if 
and only if M is nonsingular. 
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We will not prove this statement here but just remark 
that when M is singular, the matrix U is determined up 
to a unitary transformation acting on the eigenvectors of 
M with zero eigenvalue. 



1. Proof of Theorem [T] 

Thanks to the polar decomposition (|Aip . we can 
always make the order parameter A Hermitian and 
positive-semidefinite by a suitable left unitary transfor- 
mation. We thus have A = S + iA^ where S {A) is a real 
(anti)symmetric matrix. As the next and last step we 
note that a diagonal orthogonal rotation ([2]) preserves 
(anti)symmetry, and thus transforms the matrices 5, A 
separately. The real symmetric matrix S can always be 
diagonalized by such a transformation, which brings the 
order parameter to the form ([3|). 

2. Proof of Theorem [2] 

The diagonal orthogonal transformation ([2]) preserves 
Hermiticity and the spectrum, hence Eq. (jl]) says 
L/A' A, where A' = and both A and A' 

are Hermitian and positive-semidefinite. However, Theo- 
rem [7] asserts that the Hermitian part of the polar de- 
composition is unique. Therefore, we necessarily find 
A = A' = R^AR and UA = A, as was to be proved. 



APPENDIX B: INEQUALITIES BETWEEN 
QUARTIC INVARIANTS 

In this appendix we prove Theorems [3l [H [5l and [6l 
Several proofs will be based on the well known Cauchy 
inequality which asserts that for any two complex vectors 
M, w, we have 

"^VK:P). (Bl) 



The inequality in Eq. (|Bip is saturated if and only if the 
vectors u, v are coUinear. 



Proof of Theorem [ 



2. Proof of Theorem g] 

Using the complex- vector notation, C — \(j)a ■ <t>b\'^ so 
that we obviously have C > with the equality when 
all scalar products (f>a ■ <pb are zero. Writing (fia in terms 
of its real and imaginary parts, this first of all requires 
jRe^al — |Im(/)a| and Re 0a _L Im(/)a for all rows 4'a of 
the matrix A. One of the vectors, say (f)i, can then be 
cast into a special form by a right orthogonal rotation, 
so that the order parameter becomes 




Let us assume that zi ^ (if not, we go on processing 
the second row in the same manner) . The orthogonality, 
4>i'4>2 = 0, then implies Z2 + iz^ = 0, and from 4>2'4'2 = 0, 
24 = follows. The vectors (f>i and (j)2 are therefore 
coUinear so that (pi can be made zero by a suitable left 
unitary transformation. Proceeding in the same manner, 
we next nullify 4>2 and end up with just the third row of 
A, being proportional to (1,?,0). This is equivalent to 
the A-phase order parameter. 

For the second inequality in Eq. ([7]) , recall that Hermi- 
tian matrices span a real vector space with a scalar prod- 
uct defined as {A, B) = Tt:[AB). If we set Z = A^ A, we 
have B = Tr(ZZ) = ||Z||2 and C = Tr{Z* Z) = {Z*,Z). 
Therefore, we get 



< IIZ- 



\Z\? 



2{Z,Z*) = 2{B-C). 



We thus prove C < B with the equality if and only \i Z = 
Z* . This translates into the requirement that A*jAaj 
is real, i.e., all scalar products of columns of A must 
be real. However, the order parameter may always be 
transformed by a left unitary matrix to the form 




with real and complex Zj. Applying step by step the 
above reality requirement we find that the whole matrix 
actually has to be real. 



Let us denote the (real and positive) eigenvalues of 

AAt as df. Then A = {T.^^'if and B = Y.^^t■ The 
second inequality in Eq. ^ as well as the condition for 
its saturation follow immediately. For the first inequality, 
set Ui = 6f and Vi — 1. The Cauchy inequality (jBl[l then 
gives A < 3B, as was to be proved. This inequality 
is saturated when all the eigenvalues are equal, that 
is, AA^ is proportional to the unit matrix. The order 
parameter A is then unitary up to a real scale factor 
and can be brought to the CSL form by a left unitary 
transformation. 



3. Proof of Theorem [5] 

As already mentioned in Appendix IB 2[ we can define 
a scalar product of two Hermitian matrices by the trace 
of their matrix product. We now introduce some further 
notation. We first define an orthonormal basis, Ta,a = 
0, . . . , 8, so that (Ta, Tt) = Sab- These matrices are simply 
defined as Tq — t/VS and Ta — Xa/V^, a — 1, . . . , 8, in 
terms of the standard Gell-Mann matrices. Expanding a 
given Hermitian matrix in this basis, e.g. A — UaTa, the 
scalar product becomes {A, B) = aaba- 
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As above we denote Z = A. Recalling that the ba- 
sis matrices T2 , T5 , T7 are imaginary while all others are 
real, we can see that the complex conjugation in Z* just 
changes the sign of the coordinates Z2, z^, z^. Introducing 
finally the shorthand notation, 



2 2 I 2 I 2 

W — Zn + Zr. + z. 



7i 



^1 + ^3 + -^4 + ^6 + ^8 I 



we obtain the expressions for the three invariants. 



(B2) 



(TrZf = {Z,lf = 

B 

C : 



{Z, Z) = u^ + v 
{Z*,Z)^u^ + ^ 



This implies B + C = 2{u^ + v'^), whence we immediately 
get the desired inequality ([8]) . It will be saturated if and 
only if t; = 0, i.e., if the real part of Z = A^'A will 
be proportional to the unit matrix. After substitution 
for the order parameter from ([3]), it is straightforward 
to show that this is equivalent to the condition stated in 
Theorem [5l 



4. Proof of Theorem [5] 

This inequality is most tricky because it does not hold 
for all matrices A but just for those satisfying C < ^A. 
We will need another auxiliary claim. 

Theorem 8 Every Hermitian positive- semidefinite ma- 
trix Z satisfies the following inequality, 



V > wVS — u\/2, 



(B3) 



where the non-negative quantities u, v, w are defined by 
Eq. (|B2I) . The inequality is saturated if and only if Z is 
of the e type. 

Proof. — By a suitable diagonal orthogonal rotation, Z 

we can always make the coordinates Z5 and z-j 
vanish so that Z2 = ±w. Let us without lack of generality 
assume that Z2 = w. Pick up a test vector as the eigen- 
vector of A2 with the eigenvalue —1, \t) = -^(1, —h^)^ 
The expectation value of Z is 



V3 V6 \/2 

Positive-semidefiniteness requires that this be non- 
negative, which leads to 



w 



\/3 < 



^8 



L^/2 < w + uV2. 



The inequality (|B3p is thus proved. In order for it to be 
saturated, we must have v — zg, i.e., zi = z^ = z^ = 
Zfj = 0, and zg = w\/3 — u\/2. The matrix Z then reads 

form with Ai — 



f(wV3 - 



2 + 71 
w/\/2 and A2 



ma/2) As, which has the e 
— u\/3 — 



We now get to the proof of Theorem [6l For C < ^A 
we find using Theorem [5l 

2 2 2 4 

B-C>-A-2C>-A--A=-A, 



whence -^u = ^VA < \/B - C 



w V2. This means 
that the right-hand side of the inequality (|B3p is non- 
negative, and the inequality can be equivalently squared. 
This yields 



C = u"^ + -w'^ > 2w2 - 2VGu'i 



2V6UW -f 2u^, 
f 3m2 = (wV2-uV3)^ 



Vc > uVi - wV2 = Va- Vb~c. 

The saturation condition is the same as that for the in- 
equality (|B3|1 . Since the order parameter A is Hermitian 
and positive-semidefinite without lack of generality, it is 
given by the unique square root of Z = A^^ A. It is there- 
fore e-like if and only if Z is e-likc. 



APPENDIX C: GL FREE ENERGY FROM NJL 
MODEL 

In the NJL model, the mean-field free energy of a su- 
perconductor is determined by the quasifermionic exci- 
tations above the Fermi sea, and is given in terms of the 
fermion propagator (j25p as 

Here the trace is taken in the functional sense and the 
factor i comes from doubling of the number of degrees of 
freedom in the Nambu formalism. The condensate contri- 
bution, quadratic in the order parameter, is not included 
here for it essentially only serves to adjust the GL coef- 
ficient b to zero at the critical temperature. Writing the 
inverse propagator as usual as = Pg'^ + S, where the 
self-energy E is off-diagonal in the Nambu space and con- 
tains the pairing field, expanding in powers of the order 
parameter up to fourth order, and subtracting the free 
energy of the normal phase, the GL functional becomes 



T 

m 



i Tr(I?oS)' + \ MVo^f 



The traces here lead to momentum integrals which are 
evaluated in the high-density approximation ^41i] : The 
integral over the momentum magnitude, |k|, is replaced 
with one over the energy measured with respect to the 
Fermi sea, ^k, multiplied by the density of states at the 
Fermi surface, Af = 

Assuming for a moment the generic form of the gap 
matrix, ^(x) = 4>a{x)Ta, where Tq is a set of momentum- 
independent matrices in the Dirac and fermion-species 
(flavor) space, the individual terms (referred to as the 
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"gradient", "mass", and "quartic" terms with obvious 
meaning) in the GL free energy density become 



gradient term 



7C(3)P| 



327r2^2y2 



((p.qfTr(A+rtAlpr,)) 



mass term = ^AAt^^fc {Ti^A+TlKZ^n)) ^ , 
7C(3) 

quartic term = ^^^^^^ A/"0*i?!)fc0;0rf 

X (Tr(A+rtAlpr,A+TtAlpTrf))p . 



In aU the expressions, angular brackets denote averag- 
ing over directions of the indicated momentum, and the 
hats unit vectors. Also, ya,q are the Fourier compo- 
nents of the order parameter, defined as in Eq. (jl7[) . and 
the traces are now taken in the Dirac and flavor space. 
The above expressions are valid for zero as well as finite 
fermion mass and it is understood that in the projectors 
A^p the energy tp and momentum p are replaced with 
their values on the Fermi level, i.e., /i and pfP- 

One can readily check the validity of the general re- 
sult above on a particularly simple example, namely the 
(relativistic) BCS superconductor. Here we have just a 
single order parameter, 0, and the corresponding matrix 
T = 75, which ensures positive-parity pairing. All the 
Dirac traces are then equal to two and the only nontriv- 
ial angular average is ((p • q)2)p = i, so that we recover 
the well-known GL functional of the BCS theory p^ . 



•?^BCS 



7C(3)p| , 
967rVTc 



7C(3)/iPF 
327r4T2 



1. GL functional for spin-one color 
superconductors 

The calculation of the GL coefficients may be pushed 
forward by making a particular assumption on the struc- 
ture of the gap matrix. It is instructive to divide the 
calculation in two steps and assume first the spin-one 
pairing structure $ = 4>aiQaJi, where the Qa's are yet 
unspecified matrices in the fiavor space, normalized for 
the sake of convenience by TrF(Qa'9b) ~ ^ab- The bilin- 
ear terms in the GL free energy are then straightforward 
to evaluate using the explicit form of the energy projec- 
tors. 



[7] : gradient term — 



7C(3>| 
967rVr2 



1 - 



2p| 
5At2 



V • . 



(CI) 



[7]: mass term =^(^l-||)t|0,|2. (C2) 

To calculate the Dirac trace in the quartic term, we re- 
sort to the ultrarelativistic limit. In this case, the energy 
projectors reduce the spatial 7-matrices to their trans- 
verse parts, 

and the Dirac trace simplifies to 

- TrD(7_Li7±j7_Lfc7_L0 = 2{VijVki - VikVji + VaVjk)- 
Using the identity 

{'PijVki)p = ^(6(5.y4i + ^ikSji + SuSjk), 
we finally get 

[7] : quartic term = 2407r^T2 i^^ij^ki-^SikS-ji+SSuSjk) 
X ^:.^^hjrck<l>di TTF{QiQbQlQd). (C3) 

Eqs. (|C1|1 . (|C2|) . and (|C3|) summarize the general 
Ginzburg-Landau functional for a pairing with the Dirac 
structure 7^ . In particular for the spin-one pairing consid- 
ered in this paper, the gap matrix has the form (|26p . The 
fiavor trace is then Tip {QlQbQlQd) = j{SabScd + ^adhc)- 
Putting all the pieces together, we arrive at the expres- 
sion for the free energy in Eq. ((?7|) . 

While the transverse pairing can be achieved with the 
Dirac structure 7^, setting $ = 4>aiQaJoJi leads (in the 
ultrarelativistic limit) to a purely longitudinal pairing. 
In this case, the Dirac traces become trivial and we just 
quote the final result. 



[707] ■ gradient term — 



4p| 



7C(3)p| 



967rVr2 



[707] : mass term : 



[707] ■ quartic term 



3'^ 



(C4) 



MPF 

27r2 



2p| 
3^2 



t|0aP. (C5) 



7C(3)mPf 



(SijSki+SikSji+SiiSjk) 



4807r4T2 

X €^^bJrck<l>dl ^TFiQiQbQlQd)■ (C6) 



As in the previous case, Eqs. (|C4p and (|C5|) are valid for 
arbitrary fermion mass, while the quartic term (jC6[) was 
for simplicity derived in the ultrarelativistic limit. 

In the end we would like to remark that a simple spe- 
cial case of our GL functional for spin-one color super- 
conductors is the pairing of quarks of a single color and 
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two flavors in the antisymmetric flavor-singlet channel, 
studied by Buballa et al. [l3|- In this case, the flavor 

1 ic) 

matrix is given hy Q = ^T2®V^ , where T2 is the Pauh 

fc) 

matrix in flavor space and the projector on the third 
quark color. The order parameter is a complex vector, 0, 
and one immediately finds, in the ultrarelativistic limit, 

^ 2407rVTc L ' ^' ' ^' J 37r2 



It is amusing to observe that the transverse and longitu- 
dinal cases differ in the sign of the ■ </>P term, which 
results in qualitatively different forms of the ground state. 
In the 7 case the negative sign leads to a polar-like state, 
(j) ~ (0,0,1)-^. On the other hand, the positive sign in 
the 7o7 case (actually considered in [l3]) leads to an A- 
like state, (p ^ (1, i, 0)'^, with peculiar properties such as 
the existence of a single Nambu-Goldstone boson with a 
quadratic dispersion relation — the spin wave [3, [sHi • 



Kim 

4807rVT2 



2|V 



7C(3)/ipF 
9607r4T2 
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